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Abstract 

The experimental results obtained in the last few years on kaon decays 
(K— » 2-7T and, above all, Ke4 decays) allow a reliable, model independent determi- 
nation of low energy tttt scattering in the SO wave. Using them and, eventually, 
other sets of data, it is possible to give a precise parametrization of the SO wave as 
well as to find the scattering length and effective range parameter. One can also 
perform an extrapolation to the pole of the "a resonance" [/o(600)]. We obtain 
the results 



a f ( , 0) = 0.233 ± 0.013 M"\ 
and, for the a pole, 



M a = 484 ± 17 MeV, r a /2 = 255 ± 10 MeV. 
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-EXPERIMENTAL STATUS OF THE 7T7T ISOSCALAR S WAVE AT LOW ENERGY: /o(600) POLE AND SCATTERING LENGTH- 

1. Introduction. 



The lightest scalar and isoscalar mesonic "resonance", sometimes called the a resonance, has enjoyed a 
peculiar status. The Particle Data Tables (PDTW), that refer to it as the /o(600) resonance, give for its 
pole a mass between 400 and 1200 MeV, and a half width between 300 and 500 MeV; and the same values 
for the "Breit-Wigner" mass and half width: wide ranges, indeed. In fact, in the old days the resonance was 
many times reported as non-existent. On the other hand, calculations based on chiral perturbation theory 
and dispersion relations! 2,3 ! give the following values for the complex pole corresponding to it: 1 

M a = 441 ± 17 MeV, r a /2 = 230 ± 15 MeV; Rcf. 2 (1.1) 

(the errors here are purely indicative; they are obtained from the spread of the values in the articles in ref. 2) 
and 

M a = 441± 16 MeV, F a j2 = 279±? 2 5 MeV; Ref. 3. (1.2) 

Similar results are obtained by Zhou et al.,M who also use chiral perturbation theory and dispersion relations 
and get 

M a = 470 ± 50 MeV, F a j2 = 285 ± 25 MeV . 

It would be nice to be able to compare these results with experiment, without the tremendous uncertainties 
given in the PDT. Furthermore, it is desirable to use, in the comparison with experiment, only data from 
the SO wave at low and intermediate energies without introducing extraneous information on the SO wave at 
high energies, or large amounts of data from other waves. This is what is done in the present paper, where 
the only theory we use is unitarity and analyticity, the last only what follows from general field theory. 2 

In the present note it will be shown that, indeed, it is possible to get a determination of M a , r a , 
from experimental data on the SO wave alone, and with a precision comparable to the theoretical one in 
(1.1) or (1.2). In fact, we will show that, from experiment, one has the values 

M CT = 484 ± 17 MeV, r a /2 = 255 ± 10 MeV. (1.3) 

The number for M a is stable; that for F a less so. Thus, theoretical calculations appear to be reasonably 
compatible with experiment. 

Among the low energy parameters, the scattering length, and effective range parameter, 
are particularly important. They are defined as in ref. 6, 

,V2 



r Re/ ; (/) ( S ) * a«+&[V + ...; = sin£«X\ k = y/s/4-M*, 



2M 7T k 2l + 1 Jl w fe^o 

where //^ is the partial wave of definite isospin I and angular momentum I, and 5^ stands for its corre- 
sponding phase shift. The same methods that allow us to extrapolate the experimental SO amplitude to the 
sigma pole permit an extrapolation to find them. The best values we get are 

4 0) = 0.233 ± 0.013 M" 1 , 6^ 0) = 0.285 ± 0.012 M~ 3 . (1.4) 

These values are very stable; they compare well with the results from unitarized chiral perturbation theory 
(cf. the last article in ref. 2) 

4 0) = 0.23lt°;°° 3 Af-\ 4 0) = 0.300 ± 0.010 M~ 3 (1.5) 

(with only statistical errors), or with what was found from chiral perturbation theory in ref. 7, 

4 0) = 0.220 ± 0.005 M^ 1 , = 0.276 ± 0.006 M~ 3 , (1.6) 



1 The determinations (1.1) and (1.2) are, in fact, based on different methods. We quote them together to give an 
idea on the spread of theoretical calculations that, one way or the other, use chiral perturbation theory to get the 
sigma pole. 

2 In contrast, the authors of ref. 3 also use debatable theoretical information (see ref. 5); for example, on the 
(unmeasurable) scalar form factor of the pion, or on Regge theory. While in the articles in ref. 2 some theoretical 
input is used to control the left hand cut. 
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although the central value of Oq is a bit displaced. 

Finally, we can find a parametrization that represents very accurately the SO wave up to energies 
s 1 ' 2 ~ 950 MeV; see Eq. (6.1) in the Conclusions below. 

The reason why we are able to present such an accurate determination of the , , as well as 
the location of the sigma pole from experiment is the availability, in recent years, of very precise data on the 
SO wave at low energies, based on K decays, and the use of a powerful, model independent fitting technique. 

Before embarking on the details of the calculations, however, a few words have to be said on the 
meaning of the word "resonance" in connection with the a. One can give three definitions of the location of 
an elastic, background free, resonance mass and width: the point s 1 / 2 = ^ at which the phase shift crosses 
7r/2, and the width of the corresponding Breit-Wigner parametrization, r (usually called "Breit-Wigner" 
mass and width); the mass at which the energy derivative of the phase shift is a maximum and, for the width, 
the inverse of such derivative; and the pole of the partial wave amplitude in the (unphysical) Riemann sheet, 
at M a — ir a /2. Of these, the more physical one is the second: it identifies a resonance as a metastable state, 
whose lifetime is the inverse of the width (Wigner's time delay theory). 

In the case of narrow resonances, all three definitions coincide to first order in r/2M; but the 
situation for the a is very different. The SO wave phase shift for irir scattering, S^\s), docs indeed cross 
7r/2; but it does so^ at an energy of s 1 / 2 = — 800 MeV. The energy derivative of 5^\s) is nowhere 
maximum near this point, and the partial wave amplitude does not resemble a Breit-Wigner shape. Finally, 
and as we have already advanced in Eq. (1.3), there exists a pole in the second Riemann sheet, but at very 
low energy and with very large width. These are the reasons why the classification of the a as a resonance 
is so controversial. At any rate, in the present note we will be concerned with the complex, second Riemann 
sheet pole of the partial wave amplitude, without discussing the relevance of the resonance concept for it. 

2. The resonance condition 



In order to look for the pole on the second Riemann sheet associated to the sigma resonance, we will be 
using the well known fact that such a pole, located at ^J~sZ = M a — ir a /2, corresponds exactly with a zero 
of the S- matrix partial wave (s) , 

4%) = l + 2i/r( S ) = e^ 0) , 

i ii 

at sj = M a + i-Tcr/2 in the physical Riemann sheet. Thus, we can find the location of the resonance by 
looking for the solutions of SQ°\s a ) — in the upper half of the complex plane. This zero condition may be 
written in a simpler manner as 

cot^ 0) (5 ff ) = -i. (2.1) 

Solving (2.1) requires analytical continuation from the real axis. The problem with determining the 
location of the resonance pole is that analytic continuation suffers from instability problems, due to the fact 
that the function is not known exactly on the real axis. In the old days, the low energy experimental data 
for the SO wave were hopeless: several solutions existed, and the errors were very large. This is one of the 
reasons for the widely different determinations of the tr pole reported in the PDT.W Fortunately, and as 
already remarked, the situation has improved enormously in the recent years. We now have reliable data 
from just above threshold up to the kaon mass, mx, due to K; 4 decays and K 2w decays. ^ Moreover, and 
although still preliminary, we have even better very recent K; 4 decay data. I 9 ' 

In addition to this, we can profit from the fact that the determinations of the SO wave, although not 
very precise, are reasonably compatible among themselves in the range 810 - 960 MeV; see ref. 6, Eqs. (2.13). 
This is sufficient to make a fairly stable analytical extrapolation. 
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FlGURE 1. The mapping s —> w. 



3. The conformal mapping method: theory 



The unitarity of the S matrix imposes, in the elastic region, a constraint on partial wave amplitudes; for the 
SO wave, 

Im/f MH/fMI 2 

If we write 



s 1 /2^( s )/2fc-i 

with k = \J sjA. — M% the cm. momentum, then what we know of the analyticity and unitarity properties 
of the SO partial wave amplitude implies that the function (effective range function) 

V>(s) = (2fc/s 1 / 2 )cot4 0) (s) 

is analytic in the complex plane cut from — oo to and from the point where inelasticity is important, so, to 
+oo (Fig. 1): note that the elastic cut is absent in tp(s). In the case of the SO wave, one can take so = 4mf r . 

The standard mathematical method to deal with this situation is to make a conformal mapping, 
from the variable s to the variable 

y/s- \J 4m 2 K - s 

w(s) = 



y/s + v /4m|- - s 

under which the cut plane is mapped into the unit disk (Fig. 1) and the analyticity properties of the effective 
range function in the complex plane in the variable s are strictly equivalent to uniform, absolute convergence 
of a Taylor series for the function ip in the variable w in the disk \w\ < 1. For this convergence, it necessary to 
separate off those poles of the effective range function that lie on the real axis; and it can also be convenient 
to separate the zeros of the effective range function in the same region (although this last is not necessary; see 
below). Of these we have one of each: a pole due to the so-called Adler zero of the partial wave amplitude, 
lying near the left hand cut, at s = \z§, z a ~ M r (with M n the charged pion mass); and then there is a 
zero due to the phase shift crossing it/2 for an energy s 1 / 2 = (1q ~ 800 MeV. Thus, we can, in all generality, 
write the following parametrization [we will at times simplify the notation by writing S for 6^]: 

1/2 M 2 2 _ 

cot^a) = s Zl_ ™^{5 + B lW (*) + ^(s) 2 + •••}. 
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However, if using this expression with a finite number of terms, it presents the problem that a ghost is 
generated in the partial wave amplitude. 3 Although this ghost is harmless, being weak and very near the 
left hand cut (see the discussion in Appendix A) it is, as a matter of principle, better to use a ghost-free 
expansion. For this, it is sufficient to replace the formula above by 

1/2 M 2 2,2 . 

ST 7T|3 V{m^ + * + BMe> + BMs? + ■ " }■ (!U) 

Another possibility is not to separate the zero at s = yu 2 ,, writing simply 

1/2 M 2 . 2 . 

cot S(s) ^^{j^ + B^ B M s) + B 2W[S ? + •-}. 0.2) 

In these last case the zero is generated by the combination of the B n w n and we will, generally speaking, 
need one more term in the expansion than if we used (3.1), although the number of free parameters will be 
the same. 

Now, the key point for us is that the expansions (3.1), (3.2) converge in the whole cut plane: 
therefore, they can be used as they stand to solve (2.1), and hence to find the location of the resonance. 
In particular, the expansions converge at threshold and thus can also be used to determine the values of 
the scattering length and effective range parameters. As to the evaluation of the errors, we simply vary the 
parameters within the errors obtained in the fits to data, assuming that they are uncorrelated, and thus 
calculate the errors of the related observables. We have verified that, in the case of our best fit [Eq. (4.8) 
below] the parameters are actually almost uncorrelated, so this is a valid procedure. 

It is also to be noted that, because of the nature of the singularities on the cuts, we expect the B n 
to decrease like 

|£g~(n+l)- 3 / 2 (3.3) 

for large n. This last condition is, strictly speaking, only valid when one separates off the ghost piece, i.e., 
for the B n ; but, because the ghost-producing piece Zq/M^^/s is so small, it also would hold for the B n as 
well. It should also be taken into account that the behaviour (3.3) is expected to set in earlier if we separate 
out the zero, as in (3.1); if we do not separate out the zero, the first B n are constrained to cooperate to 
build up the zero, so the asymptotic regime will set later. This is in fact what we observe in our fits. 

We can get good fits to data with one or two B n s in (3.1), or two to three in (3.2); no more are 
necessary. Generally speaking, the situation is as follows: if we include in the fit only low energy (K-decay) 
data, then a x 2 /d.o.f. smaller than unity is achieved with only two parameters. If we include a third 
parameter in the fit, there is no substantial decrease of the x 2 /d.o.f., and there appear spurious minima: 
hence this third parameter is rejected by the fit. On the other hand, if we also include in the fit high 
energy data, then one can still fit with only two parameters; but if including a third parameter, the fit 
remains meaningful (the x 2 /d.o.f. decreases appreciably) and, moreover, it is theoretically more satisfactory. 
In particular, in the case where we fit with expression (3.2), the presence of terms in w, w 2 means that 
we uncorrelate the right hand cut and the left hand cut of the partial wave, while with only two B n s we 
would have an average representation of both cuts. But even if we include high energy data, one cannot 
meaningfully fit with four parameters, as spurious minima appear. All this should be clear from the values 
of the B n 's and the x 2 /d.o.f .'s of the various fits in the text. 

Another technical point is that we will fix the Adler zero to z = M n in our analysis (the fit depends 
very little on the precise value of z , provided it is near this), and we let the B n , fi vary. 

Before turning to the actual calculations a few words have to be said on the matter of analytical 
extrapolation. In principle, if one knows an analytic function on a segment of the real line, the function 
is determined on the whole (cut) plane. In practice, however, the function is not known exactly. One has, 
therefore, to test for potential instabilities in the extrapolation procedure. There are two important sources 
of instability: first, instabilities due to small variations of the central values inside the error bars of what one 
may consider the best fit. Secondly, we have the dependence of the results on the number of parameters used 
for the fits, or on the different functions used for the fits. The way to deal with this is to try extrapolations 
with fits to various data sets (provided they are compatible within errors), and to try fits with varying 

3 We are grateful to Profs. Caprini and Leutwyler for this remark. 
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FlGURE 2. The w disk, |w| < 1. The thick lines are the regions where one has 
reliable experimental data. The sigma pole is also shown; the points labeled "Ex- 
periment" are those obtained from our fits to experiment, Eqs. (4.3) and (4.10). 
The open point labeled CCL is the Caprini, Colangelo, Leutwyler calculation,' 3 ' 
with error like the thickness of the point, and the open point labeled ref. (2) is the 
average quoted in Eq. (1.2). 

number of parameters and expressions. The final errors should be enlarged to encompass the results with 
these various fits: this is what will be done here. 

4. Results 

4.1. Fits with only low energy (K decay) data 

We first fit with only data from kaon decays, generically the more reliable ones. In what follows, "Old K 
decay data" refers to the data from K14 and K 27r decays of ref. 8 (for the exact value used from the K 27r 
analysis, see Appendix B here). "NA48/2" denotes the K/4 decay data of ref. 9. For these last data statistical 
and systematic errors are taken into account; they are added quadratically. For the K/4 data of Pislak et al., 
ref. 7, we have added statistical and systematic errors in quadrature, with the exception of the higher energy 
point, that we will discuss in Appendix B. Let us note that all "Old K decay data", as well as "NA48/2" 
data, lie well inside the convergence region of the conformal expansion, as can be seen in Fig. 2. It should 
also be noted that K/ 4 decay data only give the difference between the SO and P wave phase shifts. However, 
this is not important as the P wave phase shift can be determined with great precision from the form factor 
of the pionJ 10 ! Likewise, K 27T decays give the difference between SO and S2 phases; but this last is small on 
the kaon mass, and reasonably well known there. I 11 ' 
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We start with only one parameter, Bo, plus the location of the zero, [io, m (3-1), or two parameters 
Bo, Bi in (3.2), if not separating the zero. 4 We first fit Kaon decay data, including both "Old K decay data" 
and "NA48/2". We find the following results (a °^ and b °^ are systematically given in units of M v ). 
A.I) Separating the zero, i.e., with Eq. (3.1): X 2 /d.o.f. = 15.5/(22- 2); 

B a = 17.5 ± 0.4, fi = 744 ± 24 MeV 

a 0) = 0.226 ± 0.005, 6 0) = 0.287 ± 0.008; (4.1) 
M a = 493 ± 5 MeV, r a /2 = 228 ± 9 MeV . 

A.II) Not separating the zero, i.e., with Eq. (3.2): X 2 /d.o.f. = 15.8/(22 - 2); 

B =5.1 ±0.3, Si = -18.7 ±0.6 

a 0) = 0.221 ± 0.006, 6 0) = 0.296 ± 0.008; (4.2) 
M a — 458 ± 5 MeV, r a /2 = 235 ± 3 MeV . 

The phase crosses ir/2 at the energy [i — 865 ± 8 MeV. 

It is clear that the errors in (4.1), (4.2) for the observable quantities, a °\ M a and r a /2, are 
purely statistical. While the results from (4.1) and (4.2) are similar, the differences are in general larger than 
these nominal errors. One can get a reasonable error estimate by averaging (4.1) and (4.2), and enlarging the 
errors by including a systematic uncertainty which amounts to half the difference between (4.1) and (4.2). 
In this way we get what we consider the best result using only K decay data, 

a 0) = 0.224 ± 0.006 (St.) ± 0.003 (Sys.), b { 0) = 0.292 ± 0.008 (St.) ± 0.005 (Sys.); 

M a — 476 ± 5 (St.) ± 18 (Sys.) MeV, r a /2 = 234 ± 9 (St.) ± 4 (Sys.) MeV . 4 " 

These fits permit a good determination of the low energy parameters, but, as we will see later, are less 
reliable for extrapolation to the a pole; nevertheless, we still get reasonable values for M a and F a . 

An important question is the choice of parameters. We have elected to fix the location of the Adler 
zero at zq — M w ; but we could, for example in the fit separating the zero, have fixed /io = 810 MeV (which 
is the approximate value we get from fits including data at higher energies, see below) and consider, instead 
of it, zo to be the free parameter. In this case we find 

Ai =810MeV; B = 15.7 ±0.3, z = 174±10; 

a 0) = 0.222 ± 0.009, b { 0) = 0.300 ± 0.006; (4.4) 
M a = 476 ± 5 MeV, r a /2 = 243 ± 2 MeV 

with X 2 /d.o.f. = 15.6/(22 — 2): Eq. (4.4) is quite compatible with (4.1), with an almost equal x 2 /d.o.f. For 
this reason, we will continue to fix Zq = M n in the remaining fits. 



4.2. Two parameter fits including also higher energy data 



We denote by PY05-(2.13) to the set of combined data from various experiments, at energies 810 MeV < 
g 1 / 2 < 952 MeV, collected in Eqs. (2.13) in rcf. 6 (and repeated here, Appendix B, for case of reference). 
Note that, as happened with the "K decay data", these additional data also fall well inside the convergence 
region of the conformal expansion (Fig. 2). 

4 This second method is less efficient in the present case as the data region, 4M^ < s < m 2 K , is very asymmetric with 
respect to the origin w — around which we would be expanding. 



- 6 - 



-EXPERIMENTAL STATUS OF THE 7T7T ISOSCALAR S WAVE AT LOW ENERGY: /o(600) POLE AND SCATTERING LENGTH- 

It is possible to fit all the K decay data ("Old" as well as "NA48/2"), plus PY05-(2.13), with (3.1) 
and only one Bq. We find the following results: 

B.I) If separating the zero [with (3.1)], X 2 /d.o.f. = 21.5/(31 - 2) and 

B =17-2 ±0.4, no = 771 ± 19 MeV; 

4 0) = 0.228 ± 0.005, 6 0) = 0.286 ± 0.007; (4.5) 
M a = 490 ± 5 MeV, r a /2 = 237 ± 7 MeV . 

B. II) Not separating the zero at fi , i.e, using (3.2), one finds X 2 /d.o.f. = 21.3/(31 — 2) and 

B = 2.8 ± 0.3, B x = -24.0 ± 0.6; 

a 0) = 0.213 ± 0.006, 6 0) = 0.296 ± 0.008; (4.6) 
M a = 467 ± 5 MeV, r a /2 = 213 ± 3 MeV . 

These fits with only two parameters are somewhat rigid; fortunately, it is possible to include in the 
fit one parameter more, that permits more flexibility. 

4.3. Three parameter fits including also higher energy data 

When we include the PY05-(2.13) together with all (old as well as NA48/2) K decay data, it is 
possible to include a further parameter in the fit, thus making it more flexible, getting also a better x 2 /d.o.f .: 
this was not possible with K decay data alone, since there would have been too few data and including a 
superfluous parameter would have given spurious minima. Moreover, the fits will now be more realistic 
because the terms in B 1; B 2 represent (in the average) the cuts of the effective range function, and because 
now the data are more symmetrically distributed, lying on both sides of the point around which we are 
expanding, w — 0; see Fig. 2. 

We fit as in the previous subsection and find the following results. 

C. I) If we use (3.1), i.e., we separate off the zero, the x 2 /d.o.f. becomes X 2 /d.o.f. = 21.0/(31 — 3) and we 
get 

Bo =19.0 ±0.4, Bi=4.4±0.8, ^ = 781 ± 19 MeV; 

a 0) = 0.235 ± 0.008, 6 0) = 0.282 ± 0.008; , . 

M a = 496 ± 6 MeV, B a j2 = 258 ± 8 MeV . 

C.II) Not separating the zero, i.e., with (3.2): the % 2 /d.o.f. improves significantly when introducing 
the new parameter, B 2 . We find X 2 /d.o.f. = 18.7/(31 - 3) and 

B =4.3 ±0.3, Si = -26.7 ±0.6, B 2 = -14.1 ± 1.4; 

4 0) = 0.231 ± o.oo9, 4 0) = °- 288 ± °- 009 ; ( 4 - 8 ) 

M (T =474±6MeV r a /2 = 254 ± 4 MeV . 
The zero of cot<5o°^(s) now occurs at s = with [i — 801 ± 6 MeV. 

In principle, (4.8), shown in Fig. 3, is the best fit: it takes into account the most of theoretical 
information, and it has the smallest x 2 /d.o.f. However, (4.7) also takes into account the same theoretical 
information, and its x 2 /d.o.f., while larger, is quite comparable to that of (4.8); in fact, the phase shifts 
described by (4.7) and (4.8) overlap within their errors. 

Before continuing, a few words have to be said on the errors in (4.8) [or (4.7)] . The level of precision of 
these fits is such that the errors are comparable to the expected contributions of isospin symmetry breaking, 
which should therefore be considered. For the best values of the low energy parameters and location of the 
a pole, this is not important; we will compose the results in (4.7) and (4.8), adding as systematic error half 
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300 400 500 600 700 800 900 

s" 2 (MeV) 



FIGURE 3. The data points fitted, together with the phase shifts corresponding 
to our best fit, (4.8) [full line, with error given by the dark band]. Dashed lines: 
the error band of the "old" PY05 fit in rcf. 6. The dotted line is the prediction of 
ref. 7, based on Roy equations and the value (1.6) for . In the inserted box, a 
blow-up of the low energy region. 

the difference, and the final uncertainty is substantially larger than the estimated effect of isospin breaking 
corrections; for the parametrization, we will discuss this in Sect. 6. 

As we did for the fits with only K decay data, we average the results in (4.7) and (4.8), weighted 
now each with its corresponding % 2 /d.o.f., and enlarge the errors including a systematic error of half the 
difference. In this way we find what we consider our final result, for the low energy parameters and location 
of the a pole: 



4 0) = 0.233 ± 0.010 (St.) ± 0.003 (Sys.) M" 1 , b { °> = 0.285 ± 0.009 (St.) ± 0.003 (Sys.) M~ 3 ; 
M CT =484 ±6 (St.) ±11 (Sys.) MeV, r a /2 = 255 ± 8 (St.) ± 2 (Sys.) MeV. 

Eqs. (4.9) and (4.3) are very compatible, showing the stability of the fits, and of the extrapolations 
to the sigma pole, against the number of parameters used and the formulas employed. The best values for 
the complex zero, s a are shown, in the w-plane, in Fig. 2, where we have collected the results of our fits with 
two parameters to K decay data only and their average, together with our three parameters fits. The results 
of our best estimate and the more representative of our results for the a pole and low energy parameters are 
collected in the accompanying Table (a^ and in units of M n ). 
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Process 


M CT [MeV] 


r CT /2 [MeV] 


(0) 

% 


i (0) 
n 


K decay data, 2 par., (4.1) 
K decay data, 2 par., (4.2) 


493 ±5 
458 ± 5 


228 ±9 
253 ± 3 


0.226 ± 0.005 
0.221 ± 0.006 


0.287 ±0.008 
0.296 ± 0.008 


K decay data average 


476 ± 23 


234 ± 13 


0.224 ± 0.009 


0.292 ±0.013 


K decay data + PY05, 3 par., (4.7) 
K decay data + PY05, 3 par., (4.8) 


496 ±6 
474 ±6 


258 ±8 
254 ±4 


0.235 ± 0.008 
0.230 ±0.010 


0.282 ± 0.008 
0.288 ± 0.009 


Best estimate 


484 ± 17 


255 ± 10 


0.233 ±0.013 


0.285 ±0.012 



One may wonder whether our choice of data at high energies may have biased our results. To test 
this, we will consider two alternate fits, using, instead of the PY-05 collection of data, the K decay data plus 
the data sets B, C given by Grayer et alJ 12 l We choose these solutions out of the four solutions in rcf. 12 
because they are the ones that satisfy best forward dispersion relations [apart, of course, from our solution 
(4.8) here], as shown in rcf. 6. We find X 2 /d.o.f. = 55.7/(48 - 3) and 

[K decay data plus Sol. C of Grayer et alJ 12 l] : 
B =3.57 ±0.17, B 1 = -24.3 ±0.5, B 2 = -6.3±1.3; 

roi mi (4.10a) 

a [ Q> = 0.226 ± 0.008, b { „> = 0.299 ± 0.007; 

M„ = 465 ± 5 MeV, r a = 231 ± 3 MeV . 

For solution B, we have X 2 /d.o.f. = 31.4/(41 - 3) and 

[K decay data plus Sol. B of Grayer et alJ 12 l] : 
Bo =7.63 ±0.23, B l = -23.2 ± 0.6, B 2 = -23.0 ± 1.4; 

(0) (0) ( 4 - 10b ) 

a 0) = 0.251 ± 0.011, 6 0) = 0.269 ± 0.006; 

M CT = 477 ± 7 MeV, r a = 321 ± 6 MeV . 

It should be noted that the errors in (4.10a, b) are purely nominal; in fact, the errors in the Solutions B, C of 
Grayer et al. that we have used are only the statistical errors: while the very existence of several solutions to 
the same set of raw data in ref. 12 shows that systematic errors must be substantially larger (by a factor 3 or 
more; see the comments at the end of Appendix B, and Fig. 4 here). This is reflected in the incompatibility 
of the fits. However, and in spite of these problems with the errors, these new fits show the stability of the 
results for Oq°\ and, to a lesser extent, M CT ; while the value of F a is shown to be less reliably determined. 

Finally, one may also wonder about the possibility of adding more parameters to the fits. It turns out 
that a fourth parameter would be redundant, as the x 2 /d.o.f. does not improve appreciably by introducing 
it and spurious minima appear. Specifically, we find X 2 /d.o.f. = 17.5/ (31 — 4) and the central values for the 
parameters Bq — 5.8, B\ = —22.4, B2 = —26.3 and B3 = —32.5, which do not decrease at the expected rate. 

5. The scattering length with the effective range expansion 



The scattering length (but of course, not the location of the a pole) can also be obtained with the effective 
range expansion. We expand the effective range function in powers of fc 2 , 

= — ^ + ±R k 2 /M* + ■■■■ (5.1) 

only two terms are necessary. The expansion here, however, is poorly convergent when compared to the 
conformal mapping one; due to the presence of the Adler zero at s = M 2 /2, the circle of convergence only 
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FlGURE 4. The fit to the data in Solution C of Grayer et al. [Eq. (4.10), dotted 
line], together with the error band for the "old" fit in PY05 (dashed lines), and 
our best fit here [Eq. (4.8), continuous line and gray band]. The data for Solutions 
B and C of Grayer et al. are also shown. 



extends to s = -r^M 2 , i.e., to energies s 1 / 2 < 385 MeV. To remain well inside the region of convergence of 
(5.1), we only fit data at energies s 1 / 2 < 351 MeV, and one then finds (a^ and R in units of M n ) 

X 2 /d.o.f. = 4 0) = 0.225 ± 0.007t° 15 , Ro = -[1.00 ± 0.08^ 27 ]; [Old K decay data^]. 

9 A 

X 2 / d -°- f - = a o = 0.229 ± 0.005 ±0.024, Rq = -[0.93 ± 0.05±g;||]; [NA48/2 data! 9 !] 

8 8 

X 2 /d.o.f. = — — ; 4 0) = 0.229 ±0.004 ±0.021, Rq = -[0.94 ± 0.04±g;|^]; [Old K decay±NA48/2]. 

(5.2) 

The first errors here are statistical; the second are obtained changing the point where one stops fitting to 
367 MeV or to 340 MeV. Although the results are compatible with (3.7), it is clear that the effective range 
method is much less precise and also less stable than the conformal mapping one -not surprisingly, as the 
last incorporates a lot more of analyticity information. 

6. Summary and discussion 

In summary: in the present paper we have found an accurate representation of the SO wave, given by 
the parametrization (4.8), that we consider the best of our results. However, and as noted just after this 
equation, the level of precision attained requires that we consider isospin breaking effects related, for example, 
to whether we use the parametrization for 7r + 7r~ or 7r°7r° scattering. 
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It turns out that electromagnetic corrections, which are likely the largest part of isospin breaking 
corrections, are taken into account in the analysis of systematic errors in the NA48/2 experiment. This is 
the only place where they can be important, as the errors in all other experiments are much larger than the 
estimated isospin breaking effects; for example, for the P wave, where separate fits were made in ref. 10 to 
7r°7r + and 7r + 7r~ processes, the corresponding corrections were found to be at the 1% level, except for the 
kinematical corrections. We therefore take this into account by considering that, in the parametrization in 
(4.8), that we repeat here for ease of reference, 



2k s-^Zq {Mn^s > Vs+V 4:m 2 K - s (6.1a) 



B =4.3 ±0.3, B 1 = -26.7 ±0.6, B 2 = -14.1 ± 1.4; z = M n , 
we have to interpret 



k _ f V s / 4 - M li M -n = m 7T+ = 139.57 MeV for tt+tt" scattering, 
y/s/A - M%, M-k = m^o = 134.98 MeV for 7r°7r° scattering, 



(6.1b) 



according to the process for which it is used. 

The values for scattering length, effective range parameter and location of the sigma pole from (6.1) 

are 

4 0) = 0.231 ± 0.009, 4 0) = 0.288 ± 0.009; 
M a = 474 ± 6 MeV r CT /2 = 254 ± 4 MeV . 

This is perfectly compatible, both in central values and errors, with the results reported in (4.9). Because the 
numbers in (4.9) include a systematic error due to comparison of results with two fitting formulas (something 
that is of course impossible to do for the parametrization itself) , they should be considered to be the more 
reliable results, for these quantities; so we repeat them here for ease of reference: 

4 0) = 0.233 ± 0.010 (St.) ± 0.003 (Sys.) M~\ 6 0) = 0.285 ± 0.009 (St.) ± 0.003 (Sys.) M~ 3 ; 
M a = 484 ± 6 (St.) ± 11 (Sys.) MeV, r a /2 = 255 ± 8 (St.) ± 2 (Sys.) MeV. 

These results overlap, within errors, with what was found in ref. 7 using chiral perturbation theory and 
analyticity. The value of Oq is also compatible with experimental determinations' 13 ' of this quantity, using 
a method devised by Cabibbo from K — > 7r + 7r°7r° decays, or from pionium decay' 14 ' [although our results 
arc more precise and do not depend on knowledge of a 2 ^]. 

Another question is whether one can improve our determinations. The answer is no, if we only 
use experimental data on the SO wave. But it is possible to give an independent, perhaps more precise 
determination, calculating the function cot 5^\s) for complex s with dispersion relations (Roy equations). 
This necessitates input of experimental data on other waves, and input in the high energy, Regge region, so 
it is not a determination based purely on SO wave experimental data. Preliminary results indicate that one 
would find numbers for M CT and r a /2 compatible with (1.3), and with smaller errors. This procedure could 
also give improved values for the scattering length and effective range parameters. 

Appendix A: On ghosts 



As remarked in the main text, we have separated off a term Zq/M^^/s in Eqs. (3.1), (3.2) to avoid the 
appearance of a ghost, i..c., of a spurious pole in /o°^(s) located between the Adler zero and s = 0. It turns 
out that the effect of this ghost is is negligible. To show this, we repeat our best fit, Eq. (4.8), but not 
eliminating the ghost; that is to say, we replace Eq. (3.2) by 

s 1 / 2 M 2 

cot 5( a ) = ZL_ {B + B lW (a) + B 2 w(s) 2 + •••}, (A.l) 



- n - 



-R. GARCIA-MARTIN, J. R. PELAEZ AND F. J. YNDURAIN- 



i.e., without the term that removes the ghost. Then we find (4.8) replaced by 
v 2 18 8 

WithEq. (A.l): = — ; B = 4.5 ± 0.3, Bi = -26.9 ± 0.6, B 2 = -13.5 ± 1.4; 

dof 31 — 3 

a 0) = 0.231 ± 0.010; 6 0) = 0.287 ± 0.008; 
M a = 475 ± 6 MeV r a /2 = 253 ± 5 MeV; 

this corresponds to no = 801 ± 6 MeV, and is practically indistinguishable from (4.8). Removing the ghost 
is little more than an aesthetical requirement. 

Appendix B: On experimental data 

We here say a few words have about the error we have taken for the datum at highest energy of Pislak et 
alJ 8 l In this reference, the number given is 

4 0) ((381.4 MeV) 2 ) -<Ji((381.4 MeV) 2 ) = 0.285 ± 0.014 (St.) ± 0.03 (Sys.) 

Now, this result is suspicious. The error is the smallest of all those among the data of Pislak et al.,I 8 l although 
the datum is actually an average value near the edge of phase space. Moreover, as we will see in a moment, 
the central value is incompatible with other determinations. For this reason, we have, in our fits, done as in 
ref. 6 and multiplied the statistical error by a factor 1.5: so, we have taken 

4 0) ((381.4 MeV) 2 ) - <5i((381.4 MeV) 2 ) = 0.285 ± (1.5 x 0.014). 

We could have chosen, instead of this, to add systematic and statistical error linearly for this datum, i.e., to 
fit with 

4 0) ((381.4 MeV) 2 ) - 5i ((381.4 MeV) 2 ) = 0.285 ± 0.017. 
In this case, (4.8) is replaced by 

T^TT^ S = 4.53 ±0.27, B 1 = -26.9 ± 0.6, B 2 = -15.3 ± 1.4; 
dot 31 — 3 

a 0) =0.231 ±0.010, 6 0) = 0.285 ± 0.008; Mo = 804 ± 5 MeV; 
M a = 476 ± 6 MeV T CT /2 = 258 ± 5 MeV, 

that is to say, almost identical to (4.8). 

To ascertain the degree of incompatibility of the datum of Pislak et al. with the others, we have 
repeated the fits with two other possibilities: first, adding its errors in quadrature. In this case, the \ 2 
increases to 21.6, i.e., three units above what we got in (4.8). Alternatively, if we remove the datum from 
the fit, the % 2 decreases to 15.7: the datum at 381.4 MeV of Pislak et al., with its original error, carries a 
penalty of increase of the x 2 by almost six units and would certainly bias the results if included tel quel. 
This justifies out treating its error as we have done. 

Next, we explain the value we give for the datum at m 2 K . From the decays K + — > 7r°7r + , K$ — * tt tt° 
and Ks — > 7r + 7r~ one can obtain the difference ^"'(m 2 ^) — 5 < ^'\m 2 K ). Neglecting radiative corrections, and 
with the latest results from Kloe 5 one finds 

4°Vk) - 4 2 Vk) - 51.27 ± 0.82° . 
Including radiative corrections as in ref. 15, this is corrected to 

^o°^( TO ir) — ^o\ m2 K) = 57.27 ± 0.82 (exp) ± 3 (radiative) ± 1° (chiral perturbation approximations). 



We are grateful to Dr. C. Gatti for comunicating us these results before formal publication. 
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Subtracting the value of 5^\m 2 K ) obtained from fits to experiment in ref. 6, and adding the errors linearly 
(as is advisable given the uncertainties in their evaluation) we arrive at the result we have been using in our 
fits: 

5a\ m K) =48.7 ±4.9°. (B.l) 

We have verified that a very similar result is found if using the PDTW decay data, and the prescription for 
radiative corrections of refs. 16. 

Finally and for ease of reference, we repeat here the data PY05 [ref. 6, Eqs. (2.13)] used in the fits: 

4 0) (0.870 2 GcV 2 ) =91 ±9°; 4 0) (0-910 2 GcV 2 ) = 99 ± 6° ; 4° } (0.935 2 GeV 2 ) = 109 ±8°; 
4 0) (0.912 2 GeV 2 ) =103 ±8°; 4 0) (0.929 2 GcV 2 ) = 112.5 ±13°; <5<, 0) (0.952 2 GeV 2 ) = 126 ±16°; 
4 0) (0.810 2 GeV 2 ) =88 ±6°; 4 0) (0.830 2 GcV 2 ) = 92 ± 7° ; 4 0) (0.850 2 GcV 2 ) = 94 ±6°. 

(B.2) 

Note that the errors here include systematic errors, estimated by comparing different determinations; they 
are a factor 3 or more larger than the nominal, statistical errors of the different phase shift analyses (e.g., 
those in ref. 12). 
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